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We perform microscopic calculations of the inelastic cross sections for the double and triple ex- 
citation of giant resonances induced by heavy ion probes within a semicalssical coupled channels 
formalism. The channels are defined as eigenstates of a bosonic quartic Hamiltonian constructed in 
terms of collective RPA phonons. Therefore, they are superpositions of several multiphonon states, 
also with different numbers of phonons and the spectrum is anharmonic. The inclusion of (n+1) 
phonon configurations affects the states whose main component is a n-phonon one and leads to 
an appreacible lowering of their energies. We check the effects of such further anharmonicities on 
the previous published results for the cross section for the double excitation of Giant Resonances. 
We find that the only effect is a shift of the peaks towards lower energies, the double GR cross 
section being not modified by the explicity inclusion of the three-phonon channels in the dynamical 
calculations. The latters give an important contribution to the cross section in the triple GR energy 
region which however is still smaller than the experimental available data. The inclusion of four 
phonon configurations in the structure calculations does not modify the results. 

PACS numbers: 21.60De, 21.60Jz, 24.30Cz, 25.70De 



I. INTRODUCTION 

Since their discovery (about 70 years ago) the Giant 
Resonances have been considered as the best example of 
coherent motion of nuclear systems[l]. From a macro- 
scopic point of view, this collective behavior can be con- 
sidered as high frequency, harmonic vibrations of the nu- 
clear density around its equilibrium shape. If these oscil- 
lations were harmonic then higher states of equidistant 
energies should exist. The first experimental indication 
of some structures in the excitation function in heavy 
ion collisions which could be interpreted as due to the 
population of Double Giant Resonances dates back to 
1977 (see reviews in ref. [2] and references therein). Re- 
cently there has been an unambiguous evidence of the 
existence of a Double Giant Quadrupole Resonance [3]. 
More recently also the Triple Giant Quadrupole Reso- 
nance has been observed at GANIL by using, for the 
first time, the SPEG spectrometer in conjunction with 
the INDRA 47r detector [4, 5]. Multiphonon excitations 
were clearly observed in double charge exchange reactions 
using (7r+,7r^) and (7r~,7r+) reactions[6]. The study of 
the excitation of the double GDR of very heavy nuclei 



by means of relativistic Coulomb excitation have been 
investigated in experiments performed at GSI within the 
LAND collaboration [7]. These are exclusive experiments 
where projectile fragments, neutron and gamma rays 
from the excited fragments are measured. More recently, 
the same group has found a hint for a three-phonon dipole 
state by measuring differential cross section for electro- 
magnetic fission of ^ssy relativistic energy [8]. 

The theoretical studies to have a better comprehen- 
sion of the multiphonon problematic have taken vari- 
ous directions. An extension of the quasi-particle Ran- 
dom Phase Approximation (RPA) has been used in ref. 
[9] where, in addition to the mixing of one- and two- 
phonon states, also some specific three-phonon configu- 
rations have been considered as a mechanism to generate 
the damping width of the DGDR. Other approaches ex- 
ploit the so called Brink- Axel hypothesis[10]: A Giant 
Resonance can be excited on top of any nuclear excited 
state. In this approach[ll, 12], the states to which the 
one-phonon states decay are described in terms of GOE 
(Gaussian Orthogonal Ensemble). By means of the Ran- 
dom Matrix Theory the average cross section for the ex- 
citation of the double GDR is calculated as a function 
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of the spreading and damping width. The inelastic cross 
sections of the double and triple giant dipole resonances 
in a Coulomb excitation process [13] have been also cal- 
culated. 

The excitation of collective vibrational states in heavy 
ion collisions can be viewed as due to the action of the 
mean field of each collision partner on the other. A rather 
good microscopic description of such states is given by 
Random Phase Approximation. It can be introduced as 
the lowest order in a boson expansion leading to a boson 
image of the Hamiltonian which is a sum of indepen- 
dent harmonic oscillators corresponding to each collec- 
tive mode. Thus the RPA states are pure one-phonon 
and multiphonon states and their energy is the sum of 
the energies of the single phonons. When further terms 
in the boson expansion are taken into account, anhar- 
monicities arise and the eigenstates of the Hamiltonian 
are superpositions of multiphonon states, also with dif- 
ferent numbers of phonons. 

The mean field U of each nucleus is a one body oper- 
ator. When the ground state of each nucleus is approxi- 
mated by the uncorrelated Hartree-Fock (HF) one, only 
the particle-hole (ph) part of U enters. Its bosonic image 
is linear in the phonon creation and annihilation opera- 
tors. Non linear terms arise in a natural way when the 
correlations in the ground state are taken into account 
and thus the particle-particle and hole-hole parts of U 
are no more negligible. Indeed, the bosonic images of 
such terms are non-linear in the RPA phonons. 

In several previous papers[14-16], it has been shown 
that the above mentioned anharmonicities and non- 
linearities bring to a much better agreement between the- 
oretical and experimental cross sections for the excitation 
of double giant resonance states in heavy ion collisions. 

Till now only the cross sections to one- and two-phonon 
states have been calculated within the above recalled ap- 
proach. (For simplicity, we will call "n-phonon state" an 
eigenstate of the Hamiltonian whose main component is 
a n-phonon configuration). In ref. [17], however, we have 
computed the vibrational spectrum of ^°Ca and ^°*Pb by 
extending the previous calculations with the inclusion of 
three-phonon states. It was found that the two-phonon 
states are affected by this extension. In particular, the 
double Giant Resonances (GR) are appreciably pushed 
down. The lesson we learn from this is that one can- 
not perform calculations on the double GR without tak- 
ing into account the effect of the three-phonon states. 
In view of that, in the present paper we recalculate the 
cross section for the inelastic scattering ^"Ca -|- ^°Ca at 
50 MeV/A and ^^^Pb + ^ospb at 641 MeV/A in the en- 
larged space in order to check how much our previous 
results on the cross section in the double GR region are 
modified. Then we will present calculations of the in- 
elastic cross sections for the triple excitation of the GR 
where we include the three-phonon states. In these latter 
calculations we will also take into account the effects of 
the four-phonon states. 



II. APPROACHES AND FORMALISM 

In this section we shortly recall our approach and how 
anharmonicities in the spectrum and non-linearities in 
the excitation operator are treated in it. Our model 
makes use of standard semiclassical methods techniques. 
These methods are based on the assumption that nuclei 
move on classical trajectories, while the internal degrees 
of freedom are treated quantum mechanically. 



A. Multiphonon structure and anharmonicities 

Let us denote by p (h) the single particle states which 
are unoccupied (occupied) in the HF ground state of the 
nucleus and introduce the mappings[18] 

Kh + (1 - ^)T.Bi'h'K'hBph' + ... , (1) 

p'h' 

and 

ala,,^^Bl,Bl,^; 4 a,- - ^ ^^^5^^, , (2) 

h p 

where and Bph are bosonic operators 

[Bph , -Bp' h' ] = ^pp' ^hh' ■ (3) 

The second term in the right-hand side of eq. (1) is a cor- 
rection taking care of the Pauli principle. The fermionic 
Hamiltonian is then mapped onto 

Hb = [HwB^ + HnB^B + H2aB^B^) + h.c. + (4) 
{H2iB^B^B -t- H22B^B^BB -f- H 31 B^i B^ B^i B) + h.c. 

where we have dropped indices for simplicity. The term 
Hio vanishes in the HF basis. Collective phonon opera- 
tors are introduced by means of the Bogoliubov transfor- 
mation 

Qt = J2(^;hBU-Y;M- (5) 

p,h 

The X and Y coefficients can be chosen so that the part 
of the Hamiltonian which is quadratic in the B^ and B 
operators is diagonal when expressed in terms of the 
and Q ones 

Hrpa = J2 E-QlQ- (6) 

and the X and Y satisfy the RPA equations. Of course, 
the spectrum of Hupa is harmonic. The other terms of 
the bosonic Hamiltonian (5) introduce anharmonicities 
since they mix multiphonon states among themselves. In 
our model we neglect the H31 term because it is smaller 
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than the others, as it has also been checked in an ex- 
tended Lipkin model[19]. For the remaining terms we 
keep only 

H^iB^B^B + h.c. ^ HsiQ^g^Q + h.c. , (7) 

and 

H^^B^B^BB ^ HsaQ^Q+QQ , (8) 

because the others are smaller by a factor Y/ X or powers 
of it. Therefore our bosonic Hamiltonian becomes 

+ ^22QlQlQ-^iQ'^i ■ (9) 

The eigenstates and eigenvalues of Hg are then found by 
diagonalizing it in the space of the states containing up to 
a certain number TVpho of phonons. The EII22 term mixes 
multiphonon states with the same number of phonons, 
while the H21 mixes states having number of phonons 
differing by one. In ref.[17] we have shown the results 
for ^°Ca and ^°^Pb obtained in the space with Npho=5 
and we have compared them with those of ref.[15] where 
we had Npho=2. In the next section we will mention the 
main results obtained there. 



B. Non linearities in the excitation operator 

Within a semiclassical approach to nucleus-nucleus col- 
lisions the excitation of one partner (say A) is due to the 
mean field of the other (B) acting on it. Therefore the 
excitation operator is 

W{t)=Y,Wii{t)alaj (10) 

ij 

where Wij{t) = {i\UB{R{t))\j) and Ub is the mean field 
of B (including Coulomb) which depends on time through 
the relative distance between the two nuclei. The indices 
i and j denote single particle states, both occupied and 
unoccupied, in nucleus A. When one neglects the corre- 
lation present in the ground state, only the ph terms of 
W are effective and the boson image of W is linear in 
the and Q operators. Taking into account the cor- 
relations the pp and hh parts of W cannot be neglected 
and this brings to quadratic terms 

iy = VFOO + (^M/J°Qt +^|yn,QtQ^, (11) 

V vv' 

vv' 

The first term in this equation represents the interaction 
of the two colliding nuclei in their ground state. The 

W^^ part connects states differing by one phonon, the 
W^^ term couples excited states with the same number 



of phonons, while W"^^ allows coupling between states 
differing by two phonons. These new routes of excitation 
may increase the excitation probability of the multiple 
GR. 



C. The cross section 

The inelastic scattering cross section is calculated 
within a semiclassical coupled channel approach. Let 
|$c« > denote the excited state of the nucleus of which 
we want to calculate the excitation probability. These 
states are eigenstates of the Hamiltonian (9) and there- 
fore are superpositions of multiphonon states obtained 
by diagonalizing Hq in the space containing up to N^^o 
phonon configurations. The excitation probability am- 
plitudes satisfy, for each impact parameter 6, the set of 
coupled differential equations 

a' 

(12) 

which are integrated along classical trajectories with var- 
ious impact parameters 6 [15, 16]. The cross section to 
excite the state |$a > is then calculated as 

r+oo 

aa = 2TT Pc{b)T{b)bdb. (13) 

Jo 

where Pa{b) = \Aa{b,t = +oo)\'^. The integral is over 
the whole impact parameters range which is modulated 
by the transmission coefficient T{b). 

III. EXCITATION OF *°CA 

In this section we briefly recall the results we have ob- 
tained for the structure calculation for the ^°Ca[17]. As 
stated also in the introduction, in that calculations we 
were guided by the results we got from a study of an ex- 
tended Lipkin- Meshow-Glick (LMG) model[19]. In that 
paper, the original LMG model has been extended in or- 
der to include terms that play the same role than the 
anharmonic terms of our Hamiltonian (9). The Hamil- 
tonian of such extended LMG model is still exactly solv- 
able. The relevant results can be summarized as follows: 
its diagonalization in an enlarged space including up to 
three-phonon states produces results which are very close 
to the exact ones [19]. Therefore we have followed this ap- 
proach to calculate the spectrum of ^"Ca in the space of 
one-, two- and three-phonon states. 

We have used a discrete selfconsistent HF+RPA with 
a SGII interaction, including all one-phonon states with 
J<3 which exhaust at least 5% of the EWSR. For the 
nucleus ^"Ca, we have used the nine one-phonon basis 
shown in table I. We have constructed all two- and 
three-phonon configurations out of them, without energy 
cut-off, with both natural and unnatural parity. Then 
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TABLE I: RPA one-phonon basis for the nucleus Ca. 

For each state, spin, parity, energy and percentage of the 
EWSR(isovcctor for the GDR and the IVGQR and isoscalar 
for all the other states) are reported. In the last two columns 
we report the energies of the phonons after the inclusion of 
two and three-phonon states, respectively. 



State 


.r 


E harm 


EWSR 


E-zph 








(MeV) 


(%) (MeV) (MeV) 


GMRi 


0+ 


18.25 


30 


18.36 


18.30 


GMRa 


0+ 


22.47 


54 


22.00 


21.78 


GDRi 


1" 


17.78 


56 


17.35 


17.29 


GDR2 


1" 


22.03 


10 


21.64 


21.59 


ISGQR 


2+ 


16.91 


85 


16.51 


16.44 


IVGQR 


2+ 


29.59 


26 


29.09 


29.00 


3" 


3" 


4.94 


14 


4.47 


4.40 


LEOR 


3- 


9.71 


5 


9.33 


9.28 


HEOR 


3- 


31.33 


25 


30.80 


30.89 



the Hamiltonian (9) has been diagonahzed in the space 
spanned by such states. The eigenstates arc mixed states 
whose components are of one-, two- and three-phonon 
kind 

+ E <...3l^l^2Z^3> • (14) 

Vl 1^2 1^3 

The inclusion of the three-phonon states changes the 
energies of the phonon basis of a few hundred of KeV as 
it is shown in table I. The main result of the calculation 
is that the spectrum of the two-phonon states is strongly 
modified by their coupling to the three-phonon ones. The 
diagonalization in the three-phonon space produces very 
large shifts in the energies, more than one MeV (for ''°Ca) 
in almost all the cases and always downwards. In view 
of these results, the inelastic scattering cross section to 
two-phonon states has to be recalculated and compared 
with that obtained without the inclusion of three-phonon 
configurations[16] . 

A. Cross section at 50 MeV/ A 

The semiclassical calculations for the reaction 
■*°Ca-|-'*°Ca at 50 MeV/A have been performed within 
the framework described above. In this section we discuss 
the inelastic cross section calculations done by including 
up to the three-phonon states. 

In the case we are interested in, the nuclear con- 
tribution is important. The form factors have been 
calculated [20] by employing a double folding procedure 
with the transition densities calculated within the RPA. 



Furthermore, it has been introduced an optical potential 
in order to avoid the uncertainty on the integration over 
the impact parameters. Since the optical potential takes 
into account the absorption due to all channels, we have 
introduced a procedure in order to avoid double count- 
ing the effects of the channels explicitly included in our 
calculations[16]. 

The number of two- and three-phonon states con- 
structed out of the one-phonon basis given in table I is 
more than one thousand. Considering that the ampli- 
tudes are complex quantities and that we have an equa- 
tion for each angular momentum and its projection, the 
number of time dependent coupled equations to solve 
amounts to about ten thousand. We have then to re- 
duce their number in order to render the calculation fea- 
sible. We took into account only the natural parity states 
and furthermore we have considered for the calculations 
only states with an excitation energy below 60 MeV. This 
cut off in the excitation energies guarantees that we take 
into account almost all the two-phonon states and a great 
number of the three-phonon ones. Furthermore, we took 
into account, for each state, only the components whose 
value is larger than 0.03. This choice guaranties still 
a very good normalization and reduces appreciably the 
computation time. 

The contribution of the three-phonon states to the in- 
elastic cross section for ^°Ca-|-'*°Ca at 50 MeV/A can 
be appreciated in figure 1 where the results of the cal- 
culations in the larger space (solid line) are compared 
with the two-phonon ones (dot-dashed line) [16]. Actu- 
ally, we calculated the inelastic cross section for each in- 
dividual state (14) by solving the coupled channel equa- 
tions (12). The curves presented here are always the re- 
sult of a smoothing procedure with a Lorentzian with a 
width of r=5 MeV (7 Mev for excitation energies greater 
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FIG. 1: Comparison of the cross section for Ca+ Ca at 
50 MeV/ A computed including only 1 and 2 phonons (dot- 
dashed line) with the complete calculation, namely with 1,2 
and 3 phonons (solid line) . The dashed line correspond to the 
complete calculation in the harmonic and linear case. For the 
experimental data see ref. [16] . 
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than 30 MeV) of the theoretical cross sections to the dis- 
crete levels. It appears that both the one-phonon and 
the two-phonon strengths are a little bit influenced by 
the inclusion of the three phonon states. Although the 
three-phonon configurations appreciably affect both the 
energies and wave functions of the two-phonon states[17], 
their role in the calculation of the cross section seems 
to be small. On the contrary, in the three-phonon re- 
gion the increase of the cross section is appreciable and 
one can see that the inclusion of the three phonon com- 
ponents improves the agreement with the experimental 
data although there is still some cross section missing. 
One might argue that a further enlargement of the diag- 
onalization space by including up to four-phonon states 
could reduce this discrepancy. As we will see in the next 
section this is not the case. Therefore some processes 
which are not taken into account in our approach might 
be present in this energy region. Indeed, the experimen- 
tal spectrum compared with the calculations is an exci- 
tation energy spectrum of '*''Ca in coincidence with only 
one detected proton at backward angles. This coinci- 
dence measurement makes sure that below 35 MeV no 
reaction mechanism participates to the inelastic channel. 
However, above this excitation energy, a second particle 
can be emitted in the forward direction through a dif- 
ferent reaction mechanism (such as the towing mode), 
leading to a higher cross section than for a simple exci- 
tation. This can be evidenced by means of velocity plots 
[thesis Mumu] that clearly show an asymmetry. Even- 
though one proton is emitted backward another one can 
be emitted forward contributing to an increased cross 
section around 40 MeV. 

In figure 1 it is also plotted a curve (dashed) which 
corresponds to the complete calculation in the harmonic 
and linear case. Here, once again, we would like to un- 
derline the importance of the anharmonic and non linear 
contributions. Indeed, from the figure one infers that 
their presence increases the cross section along the whole 
energy range of the calculation, expecially in the double 
and triple GR energy regions. The increase of the cross 
section in the triple GR energy region could be thought 
to be due only to the presence of many more states in 
that region when three-phonon states are included. The 
comparison with the results in the harmonic and linear 
limit, where we have the same number of states, puts in 
evidence the real origin of the strong increase. Namely, 
as already stressed in ref.[15, 16, 20], the non lineari- 
ties open new routes to the excitation of multiphonon 
states while the anharmonicities allow to populate them 
through their one- and/or two-phonon components. 

In figure 2 we compare the old two-phonon calcula- 
tion of ref.[16] with the one done including only one- 
and two-phonon channels but with the energies and wave 
functions coming from the diagonalization of the Hamil- 
tonian (9) in the space of up to three-phonon states. We 
note that the different calculations produce very similar 
results. The single resonance peak is exactly the same in 
energy, intensity and shape while a little bit of difference 



can be noted in the two-phonon region as well as in the 
high energy region. The main effect, in these regions, is 
a small shift towards lower energies of the whole curve 
which is due to the relatively strong anharmonicities we 
have found when the three-phonon states are included 
in the diagonalization. In any case, the two curves are 
almost indistinguishable also in comparison with the ex- 
perimental data. 

This study shows that the calculation of the cross sec- 
tion can be considered to be converged at the n-phonon 
level when the (n-l-l)-phonon states are taken into ac- 
count in the structure calculation but are neglected in 
the dynamical excitation. In the next section we will 
apply this recipe to the three-phonon excitation i.e. we 
include the four-phonon configurations in the structure 
but not in the coupled channel calculations. 

To get a deeper insight in the role of the three phonon 
states, in figure 3 we have decomposed the inelastic cross- 
section into the one-, two- and three- phonon compo- 
nents. The one-phonon distribution is dominated by the 
low-lying states and the giant quadrupole resonance. The 
small peak at around 30 MeV corresponds to the excita- 
tion of the HEOR state. All these states are excited dom- 
inantly at low impact parameter through the nuclear in- 
teraction. The incident energy and the projectile charge 
are not large enough to induce a strong Coulomb excita- 
tion, so the GDR cross section is four times smaller than 
the GQR one. The two-phonon contribution appears to 
be rather strong, about one order of magnitude lower 
than the single phonon component. Its structure is more 
complex, the various bumps being related to the double 
low-lying state excitation, the excitation of a GQR on 
top of the low-lying mode, the double GQR and the L=5 
component of the |GQR x HEOR > state in the high 
energy tail. The latter contribution is also visible in fig- 
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FIG. 2: Comparison of the inelastic cross-section for 

*''Ca+'"'Ca at 50 McV/A between a calculation including 
only one- and two-phonon states but with the energies and 
wave functions coining from tiie diagonalization in the space 
of up to three-phonon states (sohd hne) and the old two- 
phonon calculations (dashed line). 
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ure 4. The double GQR clearly dominates the inelastic 
spectrum around 35 MeV excitation energy. 

The three phonon component appears to be important. 

It is smaller than the two-phonon strength by a factor of 
about 3. In the high energy part it becomes the domi- 
nant contribution with a structure due to the excitation 
of a low-lying mode on top of the double GQR state and 
above 50 MeV to the triple GQR phonon. By inspection 
one can infer the difference from the old calculation where 
only one- and two-phonon states were taken into account. 
In table II we show the summed cross sections in the 3 in- 
dicated regions, around the energies corresponding to the 
excitation of one, two and three GQR phonons, respec- 
tively. One can see that the cross sections corresponding 
to the pure GQR states (within parentheses) decrease by 
a factor of about 10, each time a new phonon is excited. 
For comparison, we show also, in the second row, the 
results corresponding to the calculations done with only 
one- and two-phonon states. 

Because of the complex structure of both the two- and 
three- phonon strength the total cross-section appears 
rather smooth above the double GQR bump. However, 
the large cross-section makes possible hunting for mul- 
tiphonon states using selective signals such as specific 
decays or multipolar decomposition. 

The decomposition of the inelastic cross-section into 
various multipolcs is presented in figure 4 where we show 
the most significant contributions. In order to avoid an 
unreadable figure, full of lines, we have separated it in 
two parts. In the upper one, we show the even natural 
parity multipoles contributions while in the lower part 
we plot the odd natural parity ones. The total cross sec- 
tion is plotted in both graphs. The contribution corre- 
sponding to the angular momenta L=0 and L=l present 
appreciable peaks only in correspondence of the single 
GMR and GDR, with a contribution of the double 3~ 




' 
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FIG. 3: Decomposition of the total inelastic cross-section 
(solid line, upper ouc) for '"'Ca -|- '*°Ca at 50 MeV/A into the 
one- (dashed line), two- (solid, lower one) and three-phonon 
components (dot-dashed line). 



TABLE II: Integrated cross section (in mb) for *°Ca + *°Ca 
at 50 MeV/ A in different energy bins corresponding to one-, 
two- and threc-phonon regions, respectively. In the first row 
there are shown the results for the throo-phonou complete cal- 
culations. In the second row the results corresponding to the 
calculations done including only one and two-phonon states. 
In parenthesis the values corresponding to the single, double 
and triple GQR states. 





(14-20 MeV) 


(28-38 MeV) 


(38-60 MeV) 


3-pho 


20.83 (14.71) 


5.44 (1.34) 


2.02 (0.20) 


2-pho 


22.80 (15.62) 


5.48 (1.97) 


0.66 (-) 



(L=0 component) for the monopole case. The single 3~ 
collective state dominates the low energy region. Com- 
paring this figure with the multiphonon decomposition 
one can deduce that the bump in the 3~ strength at 
high energy, around 30 MeV, i.e. in the DGR region, is 
due to the HEOR state. The shoulder at higher energy 
is a mixing of 2 and 3 phonon states. The quadrupole 
strength presents a strong peak due to the GQR. The 
higher energy structures are mainly due to the DGQR 
state except at very high energy, above 50 MeV, which is 
dominated by the 3-GQR multiplet as can be inferred by 
the presence of the same structure in the 2+, 4+ and 6+ 
components. The strong peak in the L=5 contribution is 
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FIG. 4: Decomposition of the total inelastic cross-section 
(soUd Une, in both parts) for *°Ca+*°Ca at 50 MeV/ A into 
different angular momenta contributions. The figure is di- 
vide in two for the reader convenience. In the upper part 
we show the even natural parity ones: L=0 (dotted), L=2 
(long-dashed), L=4 (dot-dashed) and L=6 (short-dashed). In 
the lower part we plot the odd natural parity angular mo- 
menta contribution: L=l (dotted), L=3 (long-dashed), L=5 
(dot-dashed) and L=7 (short-dashed). 
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TABLE III: Results of the diagonalization for some two- 
phonon states of "'"Ca. In the first column, the states are 
labelled by their main component in the eigenvector and the 

corresponding unperturbed energy, indicated in parentheses. 
In the second column, the amplitude of the main component 
Co. Then for each total angular momentum, we show the re- 
sults of the calculation in the basis up to 2 phonon states, 
the results for the basis extended to 3 phonon states and the 
results for the basis up to 4 phonon states. The last column 
contains the results in second order perturbation theory. All 
energies are given in MeV. 



Main 


CO J'' 


< 2ph < 3ph < 4ph 




component 










order 


3"(g) 3" 


-0.91 0+ 


10.96 


9.27 


8.77 


9.20 


( 9.88) 


-0.96 2+ 


10.63 


8.89 


8.43 


8.75 




-0.96 4+ 


9.85 


8.10 


7.64 


7.96 




-0.96 6+ 


10.88 


9.12 


8.67 


8.99 




-0.92 0+ 


35.27 


33.71 


33.35 


33.59 


(35.56) 


-0.96 2+ 


35.10 


33.66 


33.33 


33.59 


Di® Qi 


0.95 1" 


34.83 


33.35 


33.05 


33.24 


(34.69) 


0.96 2" 


34.56 


33.22 


32.92 


33.16 




-0.96 3- 


34.67 


33.13 


32.82 


33.02 


Qi(8> Qi 


-0.87 0+ 


33.88 


32.47 


32.03 


32.27 


(33.82) 


0.84 2+ 


33.82 


32.47 


32.01 


32.26 




0.90 4+ 


34.02 


32.61 


32.18 


32.44 


M2(g) Di 


-0.89 1" 


40.26 


38.14 


37.72 


37.65 


(40.25) 












M2(g) Qi 


-0.73 2+ 


39.62 


37.34 


36.50 


36.80 


(39.38) 












M2(g) M2 


0.67 0+ 


45.60 


42.76 


41.15 


41.18 


(44.94) 













due to the double phonon state |GQR x 3~ > while the 
one at higher energy correspond to the excitation of the 
|GQR X HEOR > state. The 7~ strength is due to the 
double GQR build on top of the low lying 3~ state. 

Since the presented strength is built from the 
monopole, dipole, quadnipolc and octupole collective 
states, the three-phonon component goes up to L=9. 
However, it appears that the angular momenta above 
L=7, which require three-phonon excitations, are a minor 
contribution to the total strength. 



B. Role of 4 phonons 

The results from the previous calculation of anhar- 
monicities evaluated in a basis including up to three 
phonon states give us a clear indication: it is not possible 
to compute the energies of three-phonon states without 
including four-phonon states in the basis. This exten- 
sion of the basis should be sufficient if we consider the 
weak influence of the three-phonon states upon the one- 
phonon ones. Moreover the results are well reproduced 



by second order perturbation theory. The introduction of 
four-phonon states in the calculation could allow to test 
the convergence of the series, and to study their effects 
upon one- and two-phonon states. Indeed, the 3-phonon 
states would undergo a strong energy shift towards low 
energies, and their influence on 2-phonon states could be 
substantially modified. We extend our basis to 4-phonon 
states to compute anharmonicities and we will look at 
the effects on the three-phonon states. 

In order to do that, we follow the same approach de- 
scribed before: The quartic Hamiltonian matrix (9) is 
diagonalized in the space of up to four-phonon states, 
thus obtaining the new mixed eigenvectors: 

+ E C",.,.3l'^1^2i^3> (15) 

V\ y-2 1^3 l^A 

In Table III we report the results of this calculation for 
some relevant two-phonon states. Looking at the new en- 
ergies, in the sixth column, we remark that the additional 
shift imparted to the 2-phonon states by the inclusion of 
the 4-phonon states is smaller than the previous case. 
Moreover the new shift is always towards lower energies. 

The characteristics of a few 3-phonon states computed 
in the new basis are presented in Table IV. Anharmonic- 
ities of 3-phonon states are still well reproduced by sec- 
ond order perturbation theory and the correction to the 
harmonic energy is still negative. The energy shift is 
due to the presence of the 4-phonon states which push 
downwards the 3-phonon ones. This can be understood 
looking at the second order correction to the energy in 
perturbative theory, where the ratio between the matrix 
elements of the residual interaction appearing at the nu- 
merator and the difference in energy at the denominator 
determine the properties of the state. The largest matrix 
elements are those coupling 3-phonon states to 4-phonon 
ones. This is especially true in the cases involving triple 
and quadruple GMR states and arises from symmetry 
properties of the phonons, obeying the Bose statistics. 
In general, similarly to the findings of ref.[17], the ma- 
trix elements connecting a n-phonon state to that formed 
by adding to it a GMR are large. The sign of the correc- 
tion to the energy comes from the denominator, i.e. the 
difference between the energy of the considered state and 
the 4-phonon states, the latters being nearly all located 
at higher energies. 

The presence of the four-phonon states in the diago- 
nalization basis generates eigenfunctions that are more 
mixed than the ones of the previuos calculations. So we 
get states having a main component of the order of 0.5 
and several others almost as large as it. Some examples 
are given in Table IV. The extreme case is the 2+ state 
at 53.37 MeV excitation energy, whose main component 
is {Di ® £'2)2 ® Qi and appears with an amplitude of 
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TABLE IV: Results of the diagonalization in the space including up to 4-phonon states. The first column contains the name 
of the main component of the eigenvector. The harmonic energy of this state is written below the name. We show also the 
value of the c coefficient of the main component (second column) of the state as well as its parity and total angular momentum 
(third column). In the fourth and fifth columns wc show the cigciicncrgics obtained when the diagonalization is done in a space 
up to three-phonon and up to four-phonon, respectively. These results have to be compared with the result of a second order 
perturbation theory calculation (sixth column). Finally, in the last two columns, there are shown other important components 
and the corresponding c coefficient. 



Main co J" < Sph < Aph 2"'* important a 

component order components 



Ml (gi Ml 


5 Ml 


-0.499 0+ 


54.48 


53.12 


50.47 




Ml (g) 


3" (gi3 




-0.42 


54.74 
















Ml (g 


3- (gOi 


0.41 
















M -i g) M -i (g> M2 


0.26 
















Ml 


Ml 




2'.Wi 


0.22 


(Di (gDa) 


2lgQl 


-0.37 


2+ 


56.17 


53.37 


53.09 


{Di ®D2)i (g 


Qi 


0.32 


56.72 














(-Di (gC 


32)2 ® 


Oi 


-0.19 
















{Di € 


51)2)2 


(g(3" 


g)Ol)2 


-0.21 
















(3" ig 


3-)2 


(g (M2 


lg)Ql)2 


0.13 
















(Qi (g 


■)Qi)o 


(g(Oi 


® 01)2 


0.36 
















(Qi (g 


)Ql)2 


(g(Oi 


g)Oi)o 


0.34 
















(Qi (g 


3<3l)4 


(g(Oi 


lgOl)2 


0.36 
















(Qi ® Q-Cj'z 


» (Oi 


® Ol)4 


0.31 
















[Qi ® Qi)2 


Cx) (Oi 


(g Ol)2 


-0.20 


Qi (gQi (g 


Qi 


-0.91 


0+ 


50.74 


47.8 


47.3 


{Qi g 


')Ql)2 


(g(3- 


g)3-)2 


0.24 


50.73 














Qi 


(gQi 


igQi € 


5M2 


-0.27 


Qi ® Qi ig) 


Qi 


-0.91 


2+ 


50.96 


48.0 


47.5 (Qig 


)Qi)o 


(g(3- 


ig3-)2 


-0.14 


50.73 














Qi 
Qi 


(gQi 
(gQi 


® Qi € 
igQi (? 


5M2 

5 Ml 


-0.27 
-0.17 


Qi ® Qi ig) 


Qi 


-0.65 


4+ 


51.02 


48.0 


47.56 


3- 


(g3- 


® 3' $ 


5O2 


-0.18 


50.73 














Qi 


g) Qi 


(g gi » A//2 


-0.19 
















Qi 


(gQi 


(gQi g 


5 Ml 


-0.13 
















(Qiig 


)Ql)4 


(g (3- 


<g)3-)2 


-0.12 
















M2 c 


2) Di S 


5 ((32 g 


5 3-)3 


0.18 
















Ml C 


gj/Ji S 


5 (Q2 g 


5 3-)3 


0.14 


Qi (giQi ig 


Qi 


-0.92 


6+ 


51.34 


48.3 


47.85 


Qi 


(gQi 


igQi g 


5M2 


-0.27 


50.73 














(Qi CE 
Qi 


)(5l)4 

(gQi 


g) (3- 
(g)Qi € 


ig3-)2 
5 Ml 


-0.20 
-0.18 



-0.37 in the wavefunction. This state has several other 
components as large as that, which are made of 3 and 4 
phonons. 

The same is true for the L=4 state whose main com- 
ponent is Qi®Qi® Qi with an amplitude of -0.65. The 
same state is much less mixed when the diagonalization is 
done in the space up to three-phonon states. In the latter 
case its main coefficient is 0.986 with only one big com- 
ponent corresponding to the state {Di®D2)2®Qi whose 
amplitude is c = 0.15. Another interesting result, already 
found in [17] and related to the strong coupling regime, is 
the existence of some states having as second large com- 
ponent a configuration which is not directly coupled to 
the main one by the residual interaction. This is a second 
order effect due to the fact that both these configurations 
have large matrix element with another one. 



In order to check the stability of the results on the in- 
elastic scattering cross sections to two- and three-phonon 
states, we have repeated the calculations by using the en- 
ergies and wavefunctions obtained by diagonalizing the 
hamiltonian in the large space but not including the 
four-phonon channels. We dot not show the results of 
this calculations because they are almost indistinguish- 
able from the prc;vioiis ones. Once again, the inclusion 
of the n-|-l phonons is very important in the structure of 
the n-phonons but it seems it does not affect strongly the 
dynamics. Therefore we can conchide that convergence 
has been reached (at least numerically) and the discrep- 
ancy between theory and experiment in the three-phonon 
region is due to the presence in the experimental data of 
some processes which are not taken into account in our 
approach. 
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FIG. 5: Comparison of the inelastic cross-section for 
208pj,_^208pj, MeV/A computed in ref. [15] includ- 

ing only one- and two-phonon states (dashed line) and the 
complete calculation going up to three-phonons (solid). 
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FIG. 6: Decomposition of the inelastic cross-section for 
208p^_^208p^ at 641 MeV/A (sohd line) into the one- 
(dashed), two- (dot-dashed) and three-phonon components 
(dotted). 



IV. EXCITATION OF ^"**P6 

Wc have also performed calculations for multiple exci- 
tation of "^^^Ph. The inelastic cross section for the system 
^°^P6 -|- ^°^P6 has been computed for an incident energy 
of 641 MeV/ A. At this energy the nuclear contribution is 
believed to be small, so only the relativistic Coulomb ex- 
citation has been taken into account in the same way as 
it is described in ref. [15]. The collective RPA basis states 
considered in the present calculation are listed in table 
V. As in the previous case we construct all the possi- 
ble two- and three-phonon states and we diagonalize the 
Hamiltonian in this space. For this case we will not take 
into account the effects of the four-phonon states. 



TABLE V: Same as table I for 2°®Pb. 



State 


r 


E harm 


EWSR 


E2ph 








(MeV) 


(%) (MeV) (MeV) 


GMRi 


0+ 


13.61 


61 


13.42 


13.48 


GMRa 


0+ 


15.02 


28 


14.78 


14.76 


GDRi 


1" 


12.43 


63 


12.30 


12.30 


GDR2 


1" 


16.66 


17 


16.61 


16.60 


2+ 


2+ 


5.54 


15 


5.18 


5.14 


ISGQR 


2+ 


11.60 


76 


11.59 


11.55 


IVGQR 


2+ 


21.81 


45 


21.69 


21.68 


3" 


3" 


3.46 


21 


3.21 


3.19 


HEOR 


3" 


21.30 


37 


21.19 


21.20 



In figure 5 we compare the complete calculation go- 
ing up to three- phonons (solid line) with the previously 
published[15] inelastic cross-section for ^°^P6-|-^°^P6 at 
641 A/[eV/A where only one- and two- phonon states were 
included. One can see that below 35 MeV the results 



are not affected by the inclusion of three-phonon states. 
Only a small reduction of the peak around 25 MeV, cor- 
responding to the double GDR, is visible. This reduction 
can be related to the feeding of the three phonon region, 
possible in the new calculations. 

In the high energy region the contribution of the three 
phonon states appears to be important. This is confirmed 
by the decomposition of the inelastic cross-section into 
the one- , two- and three-phonon components as shown in 
figure 6. At this relativistic energies and with such heavy 
charged ions the one-phonon cross-section is clearly dom- 
inated by the GDR. In the figure, the small shoulder at 
17 MeV is due to the high lying component of the GDR 
carrying a small fraction of the dipole strength. The 
one-phonon component around 22 MeV is the isovector 
quadrupole vibration. The double phonon component 
is clearly dominated by the double-GDR excitation. In 
fact from 23 to 34 MeV the inelastic cross-section ap- 
pears to be mainly due to two-phonon states. Indeed, 
the first peak in this energy region corresponds to the 
L=2 component of \GDRi x GDRi >, the second peak 
is the L=2 component of \GDRi x GDR2 > and the 
third one is the L=3 component of \GDRi x IVGQR >. 
The double GDR2 state has a small cross section and it 
cannot be appreciated in the figure. Above 35 MeV the 
three-phonon modes provide the most important contri- 
bution to the spectrum. Indeed, the main peak is due to 
the triple DGR\ while the second one corresponds to the 
\GDRx X GDRx x GDR2 > state. 

From table VI one can see that the integrated cross- 
section for the excitation of the GDR is large at such a rel- 
ativistic energy, reaching 3.5 barns. Then a factor 10 has 
to be paid each time a new phonon is excited still leaving 
some sizeable cross-section for two and three phonon ex- 
citations. In the second row we show the results for the 
calculations done with only one- and two-phonon states. 
In the three-phonon region we gain a factor 10 when we 
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TABLE VI: Integrated cross section (in mb) for Pb+'^°^ Pb 

at 641 MeV/A in different energy bins corresponding respec- 
tively to the GDR, the two-phonon and the three-phonon re- 
gions. In tlie second row the calculations done with only one- 
and two-phonon states. In parenthesis the values correspond- 
ing to the single, double and triple GDR states. 





(8-19 MeV) 


(22-35 MeV) 


(35-45 MeV) 


3-pho 


3451.2 (3078.4) 


325.6 (227.1) 


39.2 (18.7) 


2-pho 


3510.3 (3103.3) 


348.6 (245.0) 


4.0 (-) 



introduce in the calculation the threc-phonon states. 

To get a deeper insight in the excitation process, it 
is interesting to decompose the computed inelastic spec- 
trum in various multipolaritics. Let us first start with 
the dipole strength which strongly dominates at relativis- 
tic energies (see figure 7). The GDR, which is splitted 
into a main component at 12.5 McV and a smaller peak 
around 17 McV, is of course the main contributor but 
one can observe above 35 MeV a small contribution of 
the three-GDR state coupled to spin and parity 1~. The 
quadrupolc strength is more complex. Starting at low 
energy one observes the low lying collective 2+ state fol- 
lowed by the isoscalar GQR just above 10 MeV. Except 
for the small shoulder at 22 MeV coming from the isovec- 
tor GQR, the strong bump at 25 McV can be essentially 
attributed to the double GDR coupled to 2+ . This peak 
can be directly compared with the monopole strength 
which corresponds entirely to the double GDR coupled 
to 0^. Coming back to the quadrupolc response one no- 
tices that the two phonon contribution appears as strong 
as the one phonon excitation. A multipole analysis can 
thus be an interesting way to experimentally isolate the 
multiphonon contribution. Finally, the octupole response 
presents both a 3~ and HEOR components around 3 and 
24 MeV, followed by the triple GDR state coupled to 
3^ around 35 MeV. This three phonon component cor- 
responds to the structures observed in the 1~ response 
which is nothing but the low spin member of the three- 
phonon multiplet. 

V. CONCLUSION 

In this paper we present, for the first time, microscopic 
calculations of inelastic cross sections for the triple exci- 
tation of giant resonances induced by heavy ion probes. 

We use a microscopic approach based on RPA: the mix- 
ing of three-phonon states among themselves and with 
two- and one-phonon states is considered within a bo- 
son expansion approach with Pauli corrections. This is 
equivalent to introduce anharmonicities corrections to 
the standard harmonic approximations. At the same 
time we have also introduced non-linearities in the ex- 
ternal field. 

The calculations were done by solving semiclassical 
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E (MeV) 

FIG. 7: Decomposition of the inelastic cross-section for 
208p^_^208pj, Q^-^ MeV/A (solid line) into different angu- 
lar momenta, L=0 (long-dashed), L=l (short-dashed), L=2 
(dot-dashed) and L=3 (dotted). 



coupled channel equations, the channels being superpo- 
sitions of multiphonon states. In previous calculations 

we have considered only one- and two-phonon states ob- 
taining a good agreement with the experimental cross 
section. 

In this paper we extend these microscopic calculations 
by including the three-phonon states. By diagonalizing 
a quartic microscopic Hamiltonian in the space of up to 
three-phonon states one realizes that a correct descrip- 
tion of two-phonon states requires the inclusion of one 
and three-phonon components. The anharmonicity in 
most of the cases is of the order of 1 MeV. Calculations of 
the inelastic cross section for the excitation of one-, two- 
and three-phonons states have been performed in the 
framework of this model. The cross section in the DGR 
energy region is only slightly modified. Thus the previ- 
oiisly published results arc confirmed. On the contrary, 
as one could expect, the contribution in the TGR energy 
region is quite large giving a better agreement with the 
experimental data. We have also performed calculations 
in the space of up to four-phonon states. Although the in- 
clusion of the four-phonon states is very important in the 
wave functions and energies of the three-phonon states, 
giving rise to a much stronger anharmonicity, their influ- 
ence on the dynamics is very small. 

The decomposition of the inelastic cross section into 
one-, two- and three-phonon components shows the im- 
portance of the three contributions in different region of 
the excitation energy. In the case of ^°^Pb -|- ^°^Pb at 
E/A=641 the separation in energy is very clear and one 
can distinguish the three region of interest, in the '^''Ca 
-|- ■^''Ca at 50 MeV/A case the overlap is stronger. In 
both cases we get an increase in the triple phonon en- 
ergy region showing once again the importance of the 
anharmonicity in the internal hamiltonian and the non- 
linearity in the external field. 
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